Attojoule calorimetry of mesoscopic superconducting loops 
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We report the first experimental evidence of nontrivial thermal behavior of the simplest meso- 
scopic system — a superconducting loop. By measuring the specific heat C of an array of 450,000 
noninteracting aluminum loops with very high accuracy of ~20 fJ/K, we show that the loops go 
through a periodic sequence of phase transitions (with period of an integer number of magnetic flux 
quanta) as the magnetic flux threading each loop is increased. The transitions are well described by 
the Ginzburg-Landau theory and are accompanied by discontinuities of C of only several thousands 
of Boltzmann constants fcs. 
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Halfway between micro- and macroscopic worlds, 
mesoscopic systems are known to exhibit a series of 
unique phenomena which disappear on smaller as well 
as on larger scales (see Refs. Q for review): oscil- 

lations of transition temperature in thin superconduct- 
ing cylinders 0, magnetic flux quantization magne- 
toresistance oscillations 0, persistent currents 0, etc. 
Study of these phenomena is becoming increasingly im- 
portant because mesoscopic, several nanometers in size, 
elements will be the base of this century's electronics and 
are likely to revolutionize many areas of human activ- 
ity, such as, for example, medicine, biotechnology, and 
information processing 0, HI H|- Although the electric 
transport in mesoscopic systems has been mainly studied 
during the past 20 years, the thermal transport has also 
attracted some attention very recently Yet, very 

little is known about the thermodynamic and thermal be- 
havior of mesoscopic systems (i.e., behavior of their en- 
tropy, specific heat, etc.). Meanwhile, thermodynamics 
of nanostructured systems, new phase transitions intrin- 
sic for them llll , the energies needed for their heating, 
the heat released when the system changes its state, will 
certainly be important in numerous future applications 
of nanoelectronic devices @ . 

In the present Letter we report highly sensitive spe- 
cific heat measurements performed on the simplest meso- 
scopic system — a superconducting loop of size compa- 
rable to the superconducting coherence length £(T) in a 
magnetic field. We show that even this simple system 
exhibits behavior that differs substantially from that ob- 
served in macroscopic superconductors. More precisely, 
we observe multiple phase transitions between states with 
different numbers of magnetic vortices in the increasing 
magnetic field, accompanied by discontinuities of the spe- 
cific heat as small as only few thousands of kg. These 
mesoscopic phase transitions are due to the entrance 
of superconducting vortices into the sample and have 
been recently observed in superconducting loops similar 
to ours by susceptibility measurements and by Hall 
magnetometry |l4l Il5l |. Similar phenomena have been 




CRTBT/NanoFab 




FIG. 1: (a). The suspended attojoule thermal sensor at 30 x 
magnification. The copper heater and the NbN thermometer 
can be distinguished in the upper and in the lower parts of the 
silicon membrane, respectively. The mesoscopic samples are 
deposited on the surface of the sensor between the heater and 
the thermometer by e-beam lithography, (b) . Scanning Elec- 
tron Microscope (SEM) image of the array of superconducting 
loops, (c) . SEM image of a single aluminum superconducting 
loop. 



demonstrated to exist in mesoscopic disks mis. 

Our sample is composed of 450 thousands nominally 
identical, noninteracting |2l| aluminum square loops (2 
/xm in size, w = 230 nm arm width, d = 40 nm thick- 
ness, separation of neighboring loops = 2 /im, total mass 
of the sample m = 80 ng) , patterned by electron beam 
lithography on a suspended sensor composed of a very 
thin (4 /im thick) silicon membrane and two integrated 
transducers: a copper heater and a niobium nitride ther- 
mometer (see Fig. [TJ. The setup is cooled below the 
critical temperature T c of the superconducting transition 
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by a 3 He cryostat and then its specific heat is measured 
by ac calorimetry. The technique of ac calorimetry con- 
sists in supplying ac power to the heater, thus inducing 
temperature oscillations of the thermally isolated mem- 
brane and thermometer [l9j . For the operating frequency 
in the middle of the 'adiabatic plateau' (in our case, 
the frequency of temperature modulation is / ~ 250 Hz), 
the temperature of the system 'sensor + sample' follows 
variations of the supplied power in a quasi-adiabatic man- 
ner, allowing measurements of the specific heat C with 
a resolution of SC/C > 5 x 1CP 5 for signal integration 
times of the order of one minute. This makes it possi- 
ble to measure variations of the specific heat as small as 
10 fJ/K (which corresponds to several thousands of ks 
per loop), provided that the specific heat of the sensor 
(silicon membrane, heater, and thermometer) is reduced 
to about 10-100 pJ/K. Since the temperature oscillates 
with typical amplitude of few mK, our experimental ap- 
paratus can detect energy exchanges of only few aJ at 
the lowest temperature of 0.6 K used in our experiments. 

In the absence of magnetic field, a discontinuity of the 
specific heat is observed at T c = 1.2 ± 0.02 K, corre- 
sponding to the superconducting transition of the sam- 
ple. Once the sample is in the superconducting state, 
we keep the temperature constant and turn on the mag- 
netic field H directed perpendicular to the loops' planes. 
The field is then varied very slowly (a typical run from 
to 20 mT takes 10 hours), and the specific heat C of 
the sample is measured for each value of H . Close to 
T c (above T = 0.93T C ) C exhibits monotonic behav- 
ior with H without any detectable regular fine struc- 
ture, whereas below T — 0.93T C we observe oscillations 
of C(H) with a period of approximately 0.58 mT (see 
Fig. 0) . The field H ~ 20 mT destroys superconductiv- 
ity and produces a large jump AC of the specific heat. 
The periodic character of specific heat variations is even 
more evident when looking at the Fourier transform of 
the data (see Fig. Etc)) exhibiting a well-defined peak at 
1/iJ ~ 1.72 mT -1 . The observed period of oscillations 
of 0.58 mT corresponds to the magnetic flux quantum 
$o = 2.07 x 10~ 15 Wb through a square of 1.89 fj,m side, 
which is roughly the size of our loops. At the same time, 
the amplitude of oscillations ~ 50 fJ/K corresponds to 
specific heat variations of ~ 0.1 aJ/K (~ 7500fcs) per 
loop. A similar, periodic with the external magnetic field 
H, behavior of the magnetization M of superconduc ting 
loops and disks has been recently reported in Refs. 
and ^||. From these measurements, oscillations of C 
versus H can be anticipated by using the Maxwell rela- 
tion: DC/dH oc d 2 M/dT 2 , suggesting that periodicity of 
C with H is a signature of entrance of vortices into the 
mesoscopic sample. 

A detailed description of the oscillations of C with 
magnetic field can be obtained in the framework of the 
Ginzburg-Landau (GL) theory of superconductivity [20| . 
We consider a circular loop of the same average perime- 
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FIG. 2: Specific heat data at T = 1.05 K (a) and T = 1.00 K 
(b) . The right vertical axis shows the specific heat normalized 
by its jump AG at the critical field H ~ 20 mT where the 
transition to the normal state takes place (only a part of the 
jump is shown in the figure). For these measurements, the 
current through the heater is 4 /iA and the current through 
the thermometer is 0.1 /iA. The amplitude of temperature os- 
cillations measured at the thermometer is about 5 mK. When 
the magnetic field is varied, clearly seen oscillations of C ap- 
pear with periodicity of 0.58 mT. This corresponds to the 
quantum of magnetic flux $0 threading a square of 1.89 /im 
side. (c). Fourier transform of the signal at T — 1.00 K. The 
peak at 1/H = 1.72 mT -1 corresponds to the periodicity of 
0.58 mT seen in the data. 



ter as the actual square loop and assume the GL coher- 
ence length at zero temperature £(0) ~ 0.15 //m, cor- 
responding to the mean free path £ ~ 20 nm measured 
independently. At a given magnetic flux <E> threading 
the loop, the superconducting order parameter takes the 
form tp(r) = f n (p) exp(inc/>) , where the vorticity n is a 
number of 'giant' magnetic vortices in the loop and we 
use cylindrical coordinates with the z axis perpendicular 
to the plane of the loop. The specific heat of the loop 
C = -T(d 2 F[iP(r)]/dT 2 ), where F[tp(r)} is the GL free- 
energy functional, is shown in Fig. IHIa) by dashed lines 
for different values of vorticity n and for T — 0.95T C . 
As follows from the figure, several values of vorticity are 
possible at a given value of $. In the thermodynamic 
equilibrium, the transitions between the states with vor- 
ticities n and n+1 (so-called 'giant vortex states') occur 
at <I> = (n + l/2)$o, which minimizes the free energy 
and makes the specific heat to follow the lowest of curves 
corresponding to different n in Fig. [21 This results in 
oscillations of C($) with a period of $0: m agreement 
with the experiment. The amplitude of oscillations de- 
creases when T approaches T c , and hence the oscillations 
are likely to be masked by noise at T close to T c , which 



3 




§_ o.5 m:^^fp^i^W^^ 



0.0, 



10 20 30 40 

Magnetic flux 4>/O 



FIG. 3: Specific heat of a superconducting loop computed 
from the numerical solution of the Ginzburg-Landau equa- 
tion in the increasing magnetic field and for three different 
temperatures (solid lines). Dashed lines correspond to states 
with constant vorticities (n = 0, 1, 2, . . . from left to right). $ 
is the magnetic flux threading the loop, $o is the magnetic 
flux quantum, the curves are normalized by the discontinuity 
AC of specific heat at the superconducting-to-normal transi- 
tion that takes place at $/$o — 16.5 (a), 3>/<E> — 35.5 (b), 
and $/$ ^ 45.0 (c). 
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FIG. 4: Specific heat data at lower temperatures: T = 0.85 
K (a) and T = 0.7 K (b). The right vertical axis shows 
the specific heat normalized by its jump AC at the critical 
field where the transition to the normal state takes place. At 
T = 0.85 K, the periodicity of 2 x 0.58 mT (2$ in terms 
of magnetic flux) is observed, whereas at T = 0.7 K (and at 
lower temperatures, down to T = 0.6 K) the periodicity of 
3 x 0.58 mT (3<3>o) dominates. The insets show the Fourier 
transforms of the data for < H < 16 mT and 20 < H < 28 
mT with dominant components peaked at 1/H = 0.57 mT - 
and 1/H — 1.72 mT -1 , respectively. Note the existence of 
higher harmonics due to the triangular shape of the signal at 
low field. 



explains why we did not succeed to observe any regular 
fine structure in C at T > 0.93T C . 

It is, however, not clear if during the experiment the 
superconducting loops always remain in the thermody- 
namic equilibrium as the magnetic field is varied. Once 
the flux $ surpasses (n + 1/2) <£> , the giant vortex state 
n becomes metastable and, although a state with lower 
energy (n + 1 state) exists, the latter is separated from 
the former by an energy barrier that can take a long time 
for the system to overcome . As $ is further in- 

creased, the height of the barrier decreases and finally the 
barrier vanishes. The state n then becomes unstable and 
one more giant vortex enters into the loop. To model this 
situation, we solve the GL equation for ip(r) by increasing 
the flux <I> in small steps and using the previous-step solu- 
tion as a starting point of our relaxation-type algorithm. 
The resulting dependence of the specific heat on the flux 
exhibits discontinuities at values of $ where giant vortices 
enter into the loop. The discontinuities are separated by 
A$ = $o an d the curve C(<E>) has a characteristic asym- 
metric triangular shape (solid line in Fig. I^a)). Since 
no unambiguous signature of such a shape is observed 
in the experimental data of Fig. [21 we conclude that at 
not too low temperatures (at least, above T = 0.9 K) 
our system always remains close to the thermodynamic 
equilibrium. This is probably due to the electromagnetic 
and thermal noises which provide sufficient energy for the 



system to overcome energy barriers between states with 
different vorticities well before one of the states becomes 
unstable. 

The maximum number of giant vortices that can be 
hosted by a loop of radius R and arm width w < £(T) 
can be found to be n max ~ V3-R 2 / W £,(T), yielding 
the field H max at which oscillations of C should stop: 
Umax — VS^o/^w^T). Our experiments show that the 
loop can remain superconducting up to H slightly exceed- 
ing H m ax-, but the periodicity of C with H is lost beyond 
H m ax- We believe that this is due to the transition of 
the loop from the giant vortex state for H < H ma x to 
the multivortex state for H > H ma x, as discussed, e.g., 
in Ref. 0]- At ff > H ma x, it becomes energetically fa- 
vorable for vortices (with a size of the order of £(T)) to 
enter into the width of the arms of the loop instead of 
inside the loop itself and no periodicity of C with H is 
expected. 

To study the behavior of superconducting loops at 
lower temperatures (and hence at smaller coherence 
lengths we performed the measurements of the 

specific heat at T = 0.85 K and T — 0.7 K (see Fig.©. 
In contrast to Fig. [21 discontinuities of C are now clearly 
seen, testifying that the loops do not remain in ther- 
modynamic equilibrium and explore metastable states. 
Another difference with Fig. [3 is that the periodicity 
of the signal is now a multiple of <I>o: A$ = 2$o at 
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T = 0.85 K and A$ = 3$ at T = 0.7 K. Such 'mul- 
tiple flux jumps' or 'flux avalanches' are signatures of 
simultaneous entrances of several (2 or 3) giant vortices 
into the loop and can be considered as noncquilibrium 
phase transitions between states with n and m > n + 1 
giant vortices. Such transitions are also predicted by the 
GL theory (Fig. 0Jb,c)) and have been recently observed 
in mesoscopic systems similar to ours b y u sing the bal- 
listic Hall micro-magnetometry 0, Il5l llfij. Note that 
although the observed jumps of specific heat are rather 
sharp, their width S$ ~ 0.1$o is finite. This is due to 
slightly fluctuating geometric parameters of the loops, 
which therefore exhibit transitions between different gi- 
ant vortex states at slightly different values of the applied 
magnetic field. Mutual inductance effects are expected 
to increase 5$ . It is worthwhile to note that the dis- 
continuities SC of C following from the theory are sys- 
tematically larger than the measured ones, but both the 
theory and the experiment yield very similar results for 
SC normalized by the large discontinuity AC observed at 
the transition of the sample to the normal state (see Fig. 
^for the definitions of SC and AC). SC/AC is typically 
several percent for T close to T c , but increases to ~ 25% 
for T = 0.7 K. 

The data of Fig. 01 have been obtained for the mag- 
netic field sweep rate of 2 mT/hour, a ten times faster 
run yielded identical results. However, for a ten times 
slower field sweep, finite lifetime of metastable states 
starts to play an important role, leading to a more com- 
plicated behavior of C(H) (alternating <E>o-, 2$o _ an d 
3$o-j um Ps), which is a subject of ongoing study. The 
multiple flux jumps can be understood by studying the 
stability of the state with vorticity n with respect to an 
admixture of a state with a different vorticity m > n. 
In the increasing magnetic field, the instability first oc- 
curs for m ~ n + R/ \/2£(T), favoring typical flux jumps 
m — n = 1 for T close to T c , when £(T) is large and 
R/£(T) ~ 1, and m — n > 1 at lower temperatures, 
when £(T) decreases and R/£(T) > 1. The multiple flux 
jumps are observed only in magnetic fields below H ~ 18 
mT, whereas at larger fields and up to H max the single 
flux jumps are recovered (see the insets of Fig.^b)), in 
agreement with magnetic measurements |l 1|. Theoreti- 
cal curves also show the tendency to smaller periodicity 
(A$ = 2$ instead of A$ = 3$ in Fig. EJc)) at high 
magnetic fields. 

In conclusion, the specific heat of mesoscopic super- 
conducting loops is an oscillatory function of the exter- 
nal magnetic field and exhibits discontinuities (jumps) 
at vortex entrance. Vortices enter an isolated loop one, 
two, or three at a time, depending on the temperature at 
which the experiment is performed. Definitely, the possi- 
bility of measuring the specific heat of nanoscale objects 
with very high accuracy opens quite interesting prospects 
in the field of meso- and nanoscale thermodynamics of 



superconducting as well as normal materials. 
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